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ABSTRACT 


As compared to the commonly used constant frequency radar waveforms, the 
stepped-frequency waveform can achieve high range resolution while still retaining the 
advantages of lower instantaneous receiver bandwidth and lower analog-to-digital 
sampling rate. However, the relative radial motion between the target and the stepped- 
frequency radar will result in performance degradations, such as range error, loss in 
signal-to-noise ratio, and degraded range resolution. The solution to this problem is to 
apply velocity compensation to the received signal, which can eliminate the degradations 
due to Doppler effects. Three velocity compensation schemes for the detection of a 
moving target in clutter are designed, discussed, and compared in this thesis. Also, a 
simulation is presented to verify the concepts, and simulation results are compared and 


discussed. 
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I. INTRODUCTION 


The stepped-frequency waveform is an interpulse version of pulse compression 
techniques in high resolution radars. In the stepped-frequency waveform, the carrier 
frequencies for a burst of N transmitted pulses are monotonically increased from pulse to 
pulse by a fixed frequency increment. The stepped-frequency radar can achieve high range 
resolution with the advantages of lower instantaneous receiver bandwidth and lower 
analog-to-digital (A/D) sampling rate, both of which are the major limiting factors for 
achieving high range resolution in the conventional constant frequency radar. The high 
range resolution of the stepped-frequency radar is obtained by coherently processing a 
series of N returns from the stepped-frequency pulses. The soherent processing 1s 
performed by taking the Discrete Fourier Transform (DFT) of N pulse returns, which acts 
as a matched filter when there is no relative radial motion between the radar and the 
target. The relative radial motion during the coherent processing interval will result in 
nonlinear phase shifts in ‘is received signal. As a result, the DFT processing is not 
matched to the received signal, leading to a distorted range profile. The solution to this 
problem is to apply velocity compensation to the received signal based on the target 
dandy This thesis mainly focuses on the study of velocity compensation in the 
stepped-frequency radar. The problem of - detection of a moving target in clutter will be 
discussed, which requires clutter cancellation to eliminate the clutter and velocity 


compensation to restore the target resolution. 





Three velocity compensation schemes are jess and discussed in this thesis. In 
the first scheme, the velocity compensation is performed by multiplying by a compensation 
factor in the time domain, while in the second scheme, the velocity compensation is 
performed by circular convolution with the DFT form of the compensation factor in the 
frequency domain. Both schemes are conceptually equivalent, but have different 
implementations. Their implementations will be analyzed and compared. In the first and 
the second schemes, the target velocity is assumed to be known in advance. The third 
scheme is developed based on the first scheme with the addition of a velocity 
compensation loop which enables the target with unknown velocity to be processed; 
however, the target velocity is assumed within a defined anterval. In addition, a computer 
simulation is presented to evaluate the detection of a moving target in clutter using the 
stepped-frequency waveform by means of the three velocity compensation schemes. The 
simulation results are discussed and compared. 

In Chapter II, the principle of the stepped-frequency waveform and _ its 
corresponding radar system are introduced. The discussion of high range resolution 
processing and the analysis of related radar parameters are also included. In Chapter III, 
three velocity compensation schemes for the detection of a moving target in clutter are 
investigated. The simulation results of the three compensation schemes are also presented 
in this chapter. Conclusions from the work are provided in Chapter IV. The signal 
modeling and the computer programs used in the simulation are shown in Appendices A 


and B. 





Il. PRINCIPLES OF THE STEPPED FREQUENCY RADAR 


A. STEPPED FREQUENCY WAVEFORM 


The stepped-frequency radar transmits a burst of NV pulses, whose carrier frequencies 
are monotonically increased from pulse to pulse by a fixed frequency step size, Af as 
described in Figure 2.1. Each pulse has the same pulse width, t, and is transmitted at a 
fixed pulse repetition frequency (PRF). The carrier frequency of the mth transmitted pulse 


is given by 


f,=So + G-VAS - (2.1) 
where 


f, =nominal carrier frequency 


Af = frequency step size 


The frequency spectrum constituted by a series of N returns from a target due to the 
stepped-frequency waveform is shown in Figure 2.2. Thus, the total effective bandwidth of 


the stepped-frequency waveform is given by 


Beg =NAf (2.2) 





fo = Nominal carrier frequency 


Af = Frequency step size 


3 4 N 
fot2 Af fot3 Af fo+(N -l)A f 
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Figure 2.1 The Stepped-FrequencyWaveform. 





Figure 2.2 Total Effective Bandwidth of the Stepped-Frequency 
Waveform. 











The 3-dB instantaneous bandwith of the stepped-frequency waveform is 
approximately equal to the inverse of the pulse width, 1/t, and is much less than the total 
effective bandwidth of the waveform. It is noted that the range resolution of the stepped- 
frequency waveform radar is not the conventional range resolution, ct/2. The synthetic 
range resolution achieved by coherently processing a series of N returns from a target 
depends upon the total effective bandwidth and is given as 


oe ee 2.3 
“r= 2B NAS ey 





where 


c = the speed of light, 3x10° m/sec 


The synthetic range resolution, Ar, given by Eq.2.3 for a specific target can be 
improved either by increasing the number of transmitted pulses in a burst, NV, or by 
increasing the frequency step size, Af. 

B. STEPPED FREQUENCY RADAR SYSTEM 

The implementation of the stepped-frequency radar is similar to that of a coherent 
pulsed-Doppler radar. A block diagram of a stepped-frequency radar system is shown in 
Figure 2.3. The major changes in the stepped-frequency radar are the addition of a fast 
coherent stepped-frequency synthesizer, up-conversion and down-conversion circuitry, 
and the requirement for a wideband front end [Ref.2: p.236]. The core of the system is a 


coherent stepped-frequency synthesizer with the output frequency of f,,. The frequency 
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Figure 2.3 A Block Diagram of the Stepped-Frequency Radar. 





of the synthesizer is stepped from pulse to pulse by a fixed frequency step size, Af, and is 


written as 
Ix = (n-1)Af | (2.4) 


The stepped-frequency eae is generated by first mixing the frequency of a 
coherent oscillator, f,,,, with the output frequency of the stepped-frequency synthesizer, 
Fin Following this up-conversion, the signal is converted to the final transmitted frequency 
by mixing with the frequency of a stable local oscillator, /,,.1,. The sum frequency signal is 
then pulse modulated and amplified. Thus, the carrier frequency of each transmitted pulse 
is composed of three elements: the frequency of the coherent oscillator, Foohor the frequency 
of the stable oscillator, f,,,. and the variable stepped frequency, /,,, from the stepped- 


frequency synthesizer. The carrier frequency of the mth transmitted pulse can be written as 


Sn =Feratot Sconot Sgn = So + - WAS : (2.5) 
where 


So =F eatct Soon “Nominal carrier frequency 


On the receiver side, the return signal is first mixed with the fixed frequency of the 
stable local oscillator, f,,,, and then further down-converted by mixing with a sample of 


the stepped-frequency synthesizer. The output signal after the second mixer is the 
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intermediate-frequency (IF) signal which is then passed through an IF amplifier. The 3-dB 
bandwidth of the IF amplifier is approximately equal to the inverse of the pulse width, 1/1, 
and is centered on the frequency of the coherent oscillator, f,,,. The stepped-frequency 
synthesizer is synchronized with the transmitter so that the transmitted carrier frequency 
and the output of stepped-frequency synthesizer use the same stepped frequency, /,, 
within a pulse repetition interval (PRI). As a result, multiple-time-around target echoes 
after going through the second mixer will have frequencies varying by multiples of the 
frequency step size. Therefore, the target echo whose IF signal falls outside the passband 
of the IF amplifier (i.e., target echoes from different transmitted pulses with different Pe, 
will be filtered out. The amplified signal is then split into two channels, the in-phase 
uate I, and the quadrature channel, Q, which are mixed with the frequency of the 
coherent oscillator, £,,. in the phase detector. The in-phase signal is mixed with /,,, 
directly, and the quadrature signal is mixed with 90° shifted f_,,. The outputs of the phase 
detector are sampled and quantized in the A/D converter. It should be noted that typically, 
the sampling interval is t seconds and thus each sample forms a range cell. The samples 
from the A/D converter are stored in memory and sorted by range cells. The collected NV 


samples from a range cell in which the target is located will be used for high resolution 


processing as described in Figure 2.4. 











The collected N samples from a target 


~~ ee 


f0+(N-1)Af 





Figure 2.4 N Samples from a Target are Collected and Stored in 
Memory for High Resolution Processing. 


For high PRF waveforms, a target echo may arrive in a different PRI. The output of 
the stepped-frequency synthesizer will have switched to the succeeding sented 
frequencies before receiving and processing the previously transmitted pulses. Two 
possible solutions to this problem are to switch the synthesizer to the previous transmitted 
frequencies or to have several synthesizers to sort the multiple-time-around echoes [Ref. 


2: p.269]. 


A key advantage of the stepped-frequency radar over both the ultra-wideband radar 
and the typical linear frequency-modulated (LFM) pulse compression is the narrowband 
processing which eliminates requirements for the extremely high A/D sampling rate, a 
major constraint for conventional pulse compression techniques. The disadvantage of the 
stepped frequency radar is that it takes longer time to transmit and receive N pulses for 
high resolution processing. However, in coherent radars, target returns from N transmitted 
pulses are collected before processing; thus, the stepped-frequency radar does not pose 
longer delay than the normal coherent radar. 

C. HIGH RESOLUTION PROCESSING 


The sampled output of the phase detector for the 7th pulse can be written as 
s(n) =a, exp(Jon) (2.6) 
where 


o, = the phase of the nth pulse 


a, = the magnitude of the th pulse 


The phase component of the sampled output for the wth pulse is given by 


2Rn 
Dn = 20f ie (2.7) 


where 7 is the carrier frequency of the th pulse given by Eq.2.1, and 2R,/c is the round 
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trip time of the target return due to the mth pulse. R, is the range between the target and 


the radar and is given by 


R=R, - (a-1)vT (2.8) 
where 


R, = initial target range 


T = pulse repetition interval 


vy = relative radial velocity between target and radar 


By substituting Eq.2.1 and Eq.2.8 for f and R, into Eq.2.7, we obtain 
bn =2n(fo + (n- 1p Zo @= Dr) ) 


+ ana Ro ~1)T- an fo(n —1)T 


_ Anf oRo 
ee i, 


rr ake LE n— 


7 (n—1)T (2.9) 


The expansion of the phase component $, in above equation contains four terms. 
The first term of Eq.2.9 is the constant phase resulting from the initial target range and the 
nominal carrier frequency, and is not of any importance. The second term has an induced 


frequency shift by the amount of 2A/R,/Tc, resulting from the frequency step size during 
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the round trip time of the initial target range. The third term contains a Doppler frequency 
shift of the amount 2vf, /c due to the relative radial velocity between the radar and the 
target. The final term is a frequency dispersion which is determined by the amount of 
2v(n-1)Affc due to the relative radial velocity and the output frequency of the frequency 
synthesizer. The third and fourth terms of Eq.2.9 are the phase shifts caused by the 
relative radial velocity between the radar and the target, which will result in target 
distortion in the range profiles. 

1. Generation of HRR profile 

After target samples from N pulses are collected, the high range resolution is 
obtained by coherently processing N samples. In conventional pulsed-Doppler radars, 
coherent processing is implemented by DFT as a matched filter. Similarly, the high 
resolution range (HRR) profile in the stepped-frequency radar is obtained by taking the 


N-point DFT of the collected NW target samples, which is written as 


S(k) = p> An exp (ib n)exp (- j 2a) 


pe) = | 
é Aik) k=1,2,.,N (2:10) 





N 
=) an exp(jany, 2% 
n=] 


Note that the HRR profile maps a specific range cell in which the target is located. 
Thus, the absolute range of the target is not shown in the HRR profile. Only the relative 


positions of the scatterers on the target are determined. The DFT bins in the HRR profile 
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can be interpreted as the relative synthetic range cells which resolve the relative positions 
of the target scatterers. It is also noted that the same procedure for obtaining the HRR 
profile can be applied to all other range cells. 

To reduce the sidelobe effects in the HRR profile caused by the DFT process, the 
target signal is multiplied by a weighting function prior to computing DFT. In this thesis, 
Hamming windows are employed. The weighting function can suppress the range 
sidelobes of the HRR profile; however, it also causes a degradation of the resolution as 


indicated by a scaling factor & in the following expression for the range resolution. 


_ _ke_ | 
= SNAT (2.11) 


where 


k = 1.33 for Hamming window 


Note that the DFT process functions as a matched filter for the NV target returns due 
to a burst of N transmitted pulses only in case of zero relative radial velocity between the 
radar and the target. The non-zero relative radial velocity will result in a mismatch 
between the DFT and the Doppler-shifted target returns [Ref. 3: p.4]. The result of the 
mismatch in the DFT process leads to target spreading, shifting and attenuation in the 
range profile. The target spreading and attenuation are considerably more important than 
the target shifting because they degrade the desired range resolution and decrease the 
signal-to-noise ratio. The effects of the relative radial velocity and the required velocity 


compensation will be discussed in the next section. 
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2. Maximum Unambiguous Range, KR, 

The stepped-frequency radar achieves high range resolution by generating a phase 
shift by means of linear changes in the carrier frequencies. The induced phase shift due to 
the frequency change from pulse to pulse is given by the second term of Eq.2.9 as follows 

Af 2R 
Vind = ane —(n- 1)T (2.12) 
The induced frequency shift is obtained from above equation by determining the rate of the 


phase change as 
C 


Find - za aa (2.13) 


The maximum unambiguous range, R, is obtained by substituting the maximum induced 


frequency shift equal to PRF as 


2Ry 





ay 
PRF= 7 ¢ (2.14) 
Therefore, the maximum unambiguous range, R, in the HRR profile is given by 
Ry==— (2.15) 
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The fraction of the maximum unambiguous range window R, occupied by an original 


range cell AR is given by 


AR _ cuz _ Ag (2.16) 





D. EFFECTS OF THE RELATIVE RADIAL VELOCITY 

The third and fourth terms in Eq.2.9 are the additional phase shifts due to the 
relative radial velocity between the radar and the target, resulting in a mismatch in the 
DFT process. The consequences of the mismatch in the DFT are three negative effects on 
the range profile. The first effect is target spreading in aa The second effect is the 
attenuation of target amplitude due to the spreading. The amount of spreading is 
determined by a dimensionless parameter P given by Eq.2.17 [Ref.2: p.250], which 
represents the number of the synthetic range cells (or DFT bins) that the target signal | 


spreads in the range profile during the coherent processing interval. 





p=YAT_ _vN (2.17) 


The spreading of the target signal results in a loss of synthetic range resolution, for the 
target signal has broadened from a synthetic range cell with resolution Ar up to several 


synthetic range cells. 





The third effect is the shifting of the target signal in the range profile. The target 
position will be shifted by ZL synthetic range cells, where L is given by Eq.2.18 [Ref. 3: 


p.45]. 





re| feo te WNT _ Se vN (2.18) 


Be| Ber Ar Beg Ar(PRF) 


where 


yas 


1 
fe =fo+nBef 


If the number of the shifted synthetic range = L is greater than number of DFT bins, N, 
the target signal in the range profile will be circularly shifted, resulting in a "wraparound" 
[Ref. 3: p.41]. The phenomena of the target spreading, shifting and attenuation in the 
range profile caused by the relative radial velocity are demonstrated in Figure 2.5. 

To overcome these drawbacks caused by the relative radial velocity, the received 
signal can be compensated in phase, which eliminates the effects of the relative radial 
velocity. This phase compensation depends upon the target velocity and is normally 
termed as velocity compensation factor. The velocity compensation factor will consist of 


the third and fourth terms in Eq.2.9 with 180 degrees phase rotation as given by 


comp(n) = exp (AE v(fo +(n—lAf(n—-DT ) (2.19) 
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Phenomenons of target spreading, shifting, attenuation 


100 m/sec 


200 m/sec 
| 300 m/sec 
. 1000 m/sec 


Relative range(meter) 


Magnitude(dB) 





Figure 2.5 Phenomena of Target Spreading, Shifting, and 
Attenuation for a Single Target Due to Differernt Relative Radial 
Velocities in the Range Profile. The Target Range is Unchanged. 
SNR = 13 dB, PRF =20 KHz, f, = 1GHz, N= 512, Af= 1 MHz. 
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With the application of the velocity compensation factor to the received radar signal, the 


compensated signal becomes 
C a An _ _ 
(n) = s(n)exp fs Wfo + (n- 1IAA(n—-1)T ) (2.20) 


where s(7) is given by Eq.2.6. 

For a moving target in clutter, with DFT processing, the moving target shifts 
from the clutter-filled region to the clutter-free region in the HRR profile as shown in 
Figure 2.6. However, the straightforward application of the velocity compensation factor 
to the received radar signal will simply compensate the target signal, but the clutter will 
spread in the range profile due to the unnecessary compensation velocity. Thus, clutter 
cancellation should be done before applying the velocity compensation factor in a velocity 
compensation scheme. The detailed study of the velocity compensation scheme for a 
moving target detection in clutter using the stepped-frequency waveform will be discussed 
in Chapter III. 

E. SELECTION OF PULSE WIDTH, t AND FREQUENCY STEP SIZE, Af 

The design of a stepped-frequency radar system requires a thorough understanding 
of the effects of the primary waveform parameters on the system performance. The 
selection of the primary waveform parameters, such as the pulse width, t and the 
frequency step size, Af are important in the system design. The selection of pulse width, t 


depends on the maximum target extent, &. One basic requirement for good target 
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Figure 2.6 DFT Processing Makes a Moving Target Shift to the 
Clutter-Free Region in the Range Profile. The Relative Radial 
Velocity = 500 m/sec, PRF = 20 KHz, f, = 10GHz, N= 512, Af= 
1. MHz, CNR = 35 dB, SNR =11 dB. 
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delectability is that the pulse width must encompass the maximum target extent. This 


means that the pulse width, t must be greater than or equal to the round trip time for 


length F. 
2h 
TS] — (2.21) 
or E< a (2.22) 


If the pulse width cannot encompass the target extent (in range), different scatterers on a 
target will not be located within the same original range cell, resulting in a loss of SNK and 
an incomplete HRR profile. 

The requirement in choosing the frequency step size Af is that the maximum 
unambiguous range window R, must encompass the maximum target extent to avoid 


target wraparound. That ts, 


4 = 7 > 23 
R IN, bE (2.23) 
C. 
or Af < 35 (2.24) 


When dealing with moving target detection in clutter, the maximum unambiguous range 


window, R, should be much greater than the original range cell, AR in which the target 1s 


contained. That is, 
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~ = 
R, = 2Af >> AR 5 (2.25) 
or tAf << 1 (2.26) 


Thus, the unused portion of the HRR profile is a clutter-free region which can be used to 
detect the shifted target signal. In this situation, the original range cell, AR will comprise 
only a small fraction of the maximum unambiguous range window, K,. 

When tAf = 1, the maximum unambiguous range window, KR, is equal to an 
original range cell, AR. In this situation, it is difficult to notice the target wraparound 
caused by slight changes in the target range. Furthermore, aliasing may occur if there is 
interference from adjacent range cells. If the original range cell encompasses the target and 
the clutter, the target will be camouflaged by the clutter [Ref. 1: p.168]. Thus, the 
condition of tAf= 1 is only effective for target detection in a uichiiiae environment. 

When taf > 1, the maximum unambiguous range window, R, is smaller than the 
original range cell, AR. As a result, the region of the original range cell not contained 
within the maximum unambiguous range window will fold over and result in aliasing. This 


situation must be avoided for mapping HRR profiles [Ref. 6: p.602-603]. 
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I. VELOCITY COMPENSATION IN THE STEPPED FREQUENCY RADAR 


Velocity compensation plays an important role in the stepped-frequency radar 
system, for the resolution in the range profile degrades when there is a non-zero relative 
radial velocity between the radar and the target. A target in presence of clutter can be 
clearly seen only when sufficient relative radial velocity makes the target shift into the 
clutter-free region in the range profile. With DFT processing, the moving target spreads in 
range and attenuates in magnitude. The target spreading factor P given by Eq.2.17 
determines the number of the synthetic range cells occupied by the moving target in the 
range profile. The effective synthetic resolution for a moving target is PAr instead of the 
theoretical resolution, Av. In addition to the loss of range resolution, the excessive target 
shift due to the relative radial velocity may result in target wraparound. The number of 
synthetic range cells that a moving target shifts in the range profile 1s equal to the target 
shifting factor, Z given by Eq.2.18. 

In order to overcome the drawbacks of target spreading and attenuation due to the 
relative radial motion between the radar and the target, the velocity compensation factor 
given by Eq.2.19 can be directly applied to the received radar signal. However, this will 
simply compensate the target signal, while the clutter spreads in the range profiles due to 
the unnecessary compensation velocity. Thus, in a velocity compensation scheme, the 
clutter cancellation process should be done prior to applying the velocity compensation 


factor. 
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Three velocity compensation schemes are designed and discussed in this chapter. In 
the first and the second schemes, the relative radial velocity is assumed to be known in 
advance. The third scheme deals with the target of unknown velocity, but the velocity is 
assumed within a defined interval. Computer simulations for the three compensation 
schemes are also presented in this chapter. 

A. VELOCITY COMPENSATION SCHEME 1 

The first compensation scheme for detection of moving targets in clutter is shown in 
Figure 3.1. In this scheme, the received return signal in the time domain is first windowed 
prior to taking the DFT. This DFT process generates the uncompensated range profile. At 
this stage, the clutter occupies M(tAf) synthetic range cells, which is canceled by 
multiplying by a rectangular gating function. The modified range profile is then 
transformed back to the time domain signal by an inverse DFT (IDFT). The velocity 
compensation factor is now applied to the clutter-canceled target signal in the time 
domain. After velocity compensation, the target signal is windowed and the DFT is 
performed once again. This final DFT generates a clutter-free and velocity-compensated 
HRR profile. 

Note that in this scheme, two weighting functions are employed to reduce the range 
sidelobes caused by the DFT process. The first weighting function is necessary, and the 


second weighting function gives a minor improvement in the range resolution. 
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Figure 3.1 Velocity Compensation Scheme 1. 


B. VELOCITY COMPENSATION SCHEME 2 
In this section, an alternative scheme is developed for velocity compensation. This 
new scheme is based on frequency convolution [Ref. 4: p.427] which is shown in Figure 


3.2 and written as 
[a@)ex@))=—IDFTI® +O) (3.1) 


where x,() and x,(m) are finite duration time domain sequences both having the same 
length N, while X,(k) and X,(k) are the DFTs of x,(#) and x,(7), having the same length N. 


Note that the convolution in the frequency domain is a circular convolution, not: the 


2 








common linear convolution. Thus, the result of frequency convolution will also be a 


sequence of length N. 


X1(k) 
X3(k) 


se el el 


X2(k) 
x2(n) 


*K : Circular Convolution 


x1(n) and x2(n) are time domain sequences having length N, while X1 (k) and X2(k) are 
the DFTs of x1(n) and x2(n) 





Figure 3.2 The Property of Frequency Convolution. 


In Figure 3.1, the product of x,(m) and x,(n) is replaced by the circular convolution 
of X,(k) and X,(k) and the succeeding IDFT, where x,(7) represents the output of the 
IDFT and x,(m) represent the velocity compensation factor. Thus, the first scheme 1s 
expanded as shown in Figure 3.3. The purpose of this substitution is to cancel out the two 
operations of IDFT and DFT in Figure 3.1, resulting in a modified scheme as shown in 


Figure 3.4. The cancellation of these operations is expected to speed up the computation. 
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Figure 3.3 The Expansion of Velocity Compensation Scheme |! 
According to Frequency Convolution. 
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Figure 3.4 Velocity Compensation Scheme 2. 


C. ANALYSIS AND COMPARISON OF SCHEME 1 AND SCHEME 2 

Both velocity compensation schemes are functionally equivalent. In structure, the 
major difference between the two schemes is that in scheme 1, the velocity compensation 
is done by multiplication in the time domain, while in scheme 2, the velocity compensation 
is done by circular convolution in the frequency domain. As a basis for comparison, we 
will use the “iihbet of complex multiplies required in each scheme as a measure of 
computational complexity. 

In scheme 1 (shown in Figure 3.1), there are two Hamming windows, three DFT 
processes (including IDFT), one clutter cancellation filter, and one complex multiplication 


with the velocity compensation factor. Assume that there are N samples for the receive 
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signal and N is a power of 2. The DFT process of scheme | is achieved by the Fast Fourier 
Transform (FFT) algorithm. Thus, there are N complex multiplies in each Hamming 
window, (N/2)log, N complex multiplies in each DFT process [Ref. 4: p.506], V complex 
multiplies in the clutter cancellation and N complex multiplies in the multiplication with the 
velocity compensation factor. There are a total of (3N/2)log, N + 4N complex multiplies 
involved in scheme 1 provided that an FFT algorithm is employed in the DFT process. For 
example, when N = 512, there will be a total of 8960 complex multiplies required in 
scheme 1. If the DFT process is not implemented by an FFT algorithm, but by the 
standard DFT algorithm, there will be NV * complex multiplies in each DFT process [Ref. 4: 
p.426]. Thus, there will be a total of 3V 7 + 4N complex multiplies involved in scheme 1. 
When N = 512, there are a total of 788,480 complex multiplies required (approximately 88 
times the computational cost of using the FFT algorithm). Therefore, the FFT algorithm 
used in the DFT process of scheme 1 can provide much higher computational speed than 
the standard DFT algorithm. 

In scheme 2 (shown in Figure 3.4), we can find one Hamming window, two DFT 
processes, one clutter cancellation filter and one circular convolution with the velocity 
compensation factor in the frequency domain. A direct approach for the circular 
convolution process requires N ? complex multiplies [Ref. 4: p.426]. Assume that the FFT 
algorithm is used in the DFT process of scheme 2. Thus, there are a total of (V log, NV) + 
2N + N? complex multiplies. When NV = 512, there will be 267,776 complex multiplies 


required in scheme 2, which are much more than those required for the FFT algorithm in 
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scheme 1. If the standard DFT algorithm is employed ‘ scheme 2, there will be a total of 
3N ? + 2N complex multiplies. Thus, for N = 512, there are 787,456 complex multiplies 
required, which is about the same computational complexity as using the standard DFT 
algorithm in scheme 1. As expected, for both velocity compensation schemes, the FFT 
algorithm will offer less computational complexity and higher processing speed than the 
standard DFT algorithm. 

From the previous analysis, it is found that when the FFT algorithm is employed in 
the DFT process of both schemes, the circular convolution will be the major factor that 
results in the much greater computational complexity in scheme 2. Therefore, scheme 1 
mechanized with the FFT algorithm is recommended for velocity compensation in the 
stepped-frequency radar. 

D. VELOCITY COMPENSATION SCHEME 3 

The third scheme is proposed to deal with the target of unknown velocity. 
However, the target velocity is assumed within a defined interval. Scheme 3, as shown in 
Figure 3.5 is implemented by using the model of scheme 1 or scheme 2, with the addition 
of a velocity compensation loop. Scheme 1 is employed here as the base model because it 
has much higher processing speed than scheme 2. In scheme 3, the stepped-frequency 
radar will automatically try several compensation velocities to generate compensated HRR 
profiles. These compensated HRR profiles will be stored in memory and compared to find 
the one with the best resolution and signal-to-noise ratio (SNR). It should be noted that in 


scheme 3, the compensation velocities which begin at a defined initial value are 
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Figure 3.5 Velocity Compemsation Scheme 3. 
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monotonically increased by a velocity step size, Av in each loop. The selection of velocity 
step size, Av and the number of compensation loops will depend on the size of the defined 
velocity interval and the waveform parameters. 
E. SIMULATION 

Scenario: A moving target approaches a surface-based stepped-frequency radar at 
nose aspect as shown in Figure 3.6. The radar platform and the clutter are assumed 
stationary. The target with four dominant scatterers on it, is 8 m long at a distance of 
11250 m from the radar. The ranges of the scatterers on the target are 11251 m for 
scatterer 1] (RCS = 0.6 m’), 11253 m for scatterer 2 (RCS = 0.5 m’ ), 11255 m for 
scatterer 3 (RCS = 1 m’) and 11257 m for scatterer 4 (RCS = 0.7 m’). The relative radial 
velocity between the target and the radar is initially assumed to be 750m/sec, and will be 
varied in the simulation for the purpose of testing the effects of the relative radial velocity. 
The clutter-to-noise ratio (CNR) is 35 dB. The variance of the thermal noise which is 
modeled as white Gaussian noise is set to 0.0002. The signal-to-noise ratio (SNA), 13 dB 
is based on the power ratio of scatterer 3 (RCS = 1 m’ ) to the thermal noise. 


Waveform parameters used in the simulation are as follows. 


Number of pulses per burst (1): 512 
Nominal carrier frequency (/,): 10 GHz 
Frequency step size (Af): 1 MHz 
Pulse repetition frequency (PRF): 200 KHz 
Pulse width (7): | 0.1 psec 
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The fraction of R, occupied by AR: 


AR ct/2 _ _ 
Ro wy 


The number of synthetic range cells within the fraction of AR in the HRR profile: 








o x N= (tAf)x N=0.1x 512 =51.2%52 


Frequency resolution of DFT: 


1 _PRF_ 200x103Hz_ 4, 
an or oe = 390.625 Hz 





The maximum unambiguous range, R, of the HRR profile and its corresponding 
maximum induced frequency, PRF, in the DFT is described in Figure 3.7. Simulation 
programs are written in MATLAB where the DFT process is computed by using a radix-2 


FFT algorithm when the length of the sequence is a power of 2 [Ref.5: p.271]. 


PRE /N = 390.625 Hz 
<-- 





Figure 3.7 The Maximum Unambiguous Range, R, in the HRR 
Profile is Equivalent to the Maximum Induced Frequency Shift, 
PRF. 
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1. Effects of The Relative Radial Velocity 

To demonstrate effects of the relative radial velocity in the range profiles clearly, the 
relative radial velocity is varied in each uncompensated range profile. The relative radial 
velocities to be tested in this simulation are 0 m/sec, 350 m/sec, 650 m/sec, 950 m/sec and 
1250 m/sec. The simulation results are shown in Figure 3.8 through Figure 3.12. 

In Figure 3.8, the target signal is in clutter and cannot be detected since there is no 
relative radial velocity between the target and the radar. In Figure 3.9, the target signal 
from a moving target shifts from the clutter-filled region to the clutter-free region, but the 
magnitude of target signal is much smaller than that of the clutter. Thus, the target will 
possibly be treated as noise and neglected. As shown in Figure 3.10 through 3.12, the 
scatterers on the target cannot be resolved, and the phenomena of the target spreading, 
shifting, and signal attenuation in the range profiles become more prominent as the relative 
radial velocity increases. It is also noted that in Figure 3.9, the relative radial velocity _ 
350 m/sec, but the scatterers on the target can still be resolved because the use of a high 


PRF (200 KHz) reduces target spreading. 
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Figure 3.8 The Uncompensated HRR Profile for a Moving Target 
with Scatterers in Clutter. The Target are Masked by the Clutter 
Since There -is No Relative Radial Velocity between the Radar and 
the Target. 
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Figure 3.9 The Uncompensated HRR Profile for a Moving Target 
with Scatterers in Clutter. The Relative Radial Velocity = 350 
m/sec. 
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Figure 3.10 The Uncompensated HRR Profile for a Moving Target 
with Scatterers in Clutter. The Relative Radial Velocity = 650 
m/sec. 
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Figure 3.11 The Uncompensated HRR Profile for a Moving Target 
with Scatterers in Clutter. The Relative Radial Velocity = 950 
m/sec. 
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Figure 3.12 The Uncompensated HRR Profile for a Moving Target 
with Scatterers in Clutter. The Relative Radial Velocity = 1250 
m/sec. : 


2. Simulation of Scheme 1 And Scheme 2 

In this subsection, velocity compensation scheme 1 and scheme 2 designed for the 
detection of a moving target in clutter are simulated. The radar parameters given in the 
scenario are used in the simulation. The relative radial velocity between the radar and the 
target is 750 m/sec, and the PRF is 200 KHz. Several compensation velocities, such as 
150 m/sec, 350 m/sec, 550 m/sec, 750 m/sec, 1000 m/sec, and 1250 m/sec, are used to 
evaluate the performance of both compensation schemes. When the compensation velocity 
is equal to 750 m/sec, the range profile is correctly compensated. The uncompensated 
range profile is shown in Figure 3.13. The simulation results for compensated profiles are 


shown in Figure 3.14 through Figure 3.25. 
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In Figure 3.14 through Figure 3.25, the clutter has been canceled and the target 
resolution becomes finer as the compensation velocities get closer to the relative radial 
velocity, 750 m/sec. In the correctly compensated HRR profiles shown in Figure 3.17 and 
Figure 3.22, the scatterers on the target are the most clearly resolved and have the best 
SNR, for the relative radial velocity is fully compensated, and hence, the DFT process 
functions as a matched filter to the received signal. Also, the SNR (based on the third 
target scatterer to thermal noise) in the fully compensated HRR profiles is almost up to 30 
dB which is much higher than the previously given 13 dB. This increased SNR is due to 
the coherent integration of NV pulses in the DFT process. 

Note that in some of the HRR profiles without correct velocity compensation, for 
example, those shown in Figures 3.16, 3.18 , 3.21, and 3.23, the scatterers on the target 
can still be distinguished because the PRF used in the simulation has a high value of 200 
KHz which leads to lower target dispersion. It is also observed from Figures 3.18 and 3.23 
that the target signal is circularly shifted to the other edge in the range profile due to the 
excessive compensation velocity. 

In general, the first scheme and the second scheme generate very similar 
compensated HRR profiles for the same compensated velocity. Nevertheless, comparing 
the simulation results of both schemes carefully, it is found that the compensated HRR 
profiles generated by scheme 2 have higher range sidelobes and noise ripples. The lower 
range sidelobes in the compensated HRR profiles generated by scheme 1 are due to the 


existence of one more Hamming window in the scheme. 
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Figure 3.13 The Uncompensated HRR Profile. The Relative Radial 
Velocity = 750 m/sec. 
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Figure 3.14 The Compensated HRR Profile of Scheme 1. The 
Clutter is Canceled. Compensation Velocity = 150 m/sec. 
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Figure 3.15 The Compensated HRR Profile of Scheme 1. The 
Clutter is Canceled. Compensation Velocity = 350m/sec. 
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Figure 3.16 The Compensated HRR Profile of Scheme 1. The 
Clutter is Canceled. Compensation Velocity = 550m/sec. 
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Figure 3.17 The Compensated HRR Profile of Scheme 1. The 
Clutter is Canceled. Compensation Velocity = 750 m/sec. Note 
That This is a Correctly Velocity-Compensated HRR Profile. 
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Figure 3.18 The Compensated HRR Profile of Scheme 1. The 
Clutter is Canceled. Compensation Velocity = 1000 m/sec. The 
Target Scatterers are Circularly Shifted to Another Edge of the 
Range Profile Due to the Excessive Compensation Velocity. 
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Figure 3.19 The Compensated HRR Profile of Scheme 1. The 
Clutter is Canceled. Compensation Velocity = 1250 m/sec. The 
Target Scatterers are Circularly Shifted to Another End of the 
Range Profile Due to the Excessive Compensation Velocity. 
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Figure 3.20 The Compensated HRR Profile of Scheme 2. The 
Clutter is Canceled. Compensation Velocity = 150 m/sec. 
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Figure 3.21 The Compensated HRR Profile of Scheme 2. The 
Clutter is Canceled. Compensation Velocity = 350 m/sec. 
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Figure 3.22 The Compensated HRR Profile of Scheme 2. The 
Clutter is Canceled. Compensation Velocity = 550 m/sec. 
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Figure 3.23 The Compensated HRR Profile of Scheme 2. The 
Clutter is Canceled. Compensation Velocity = 750 m/sec. This is a 
Correctly Velocity-Compensated HRR Profile. 
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Figure 3.24 The Compensated HRR Profile of Scheme 2. The 
Clutter is Canceled. Compensation Velocity = 1000 m/sec. The 
Target Scatterers are Circularly Shifted to Another End of the HRR 
Profile Due to the Excessive Compensation Velocity. 
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Figure 3.25 The Compensated HRR Profile of Scheme 2. The 
Clutter is Canceled. Compensation Velocity = 1250 m/sec. The 
Target Scatterers are Circularly Shifted to Another End of the HRR 
Profile Due to the Excessive Compensation Velocity. 





3. Simulation of Scheme 3 

In the simulation of scheme 3, the target approaches the radar platform at nose 
aspect at a velocity between Mach 2 and Mach 4 (approx. 660 m/sec to 1320 m/sec). The 
given initial compensation velocity is 660 m/sec, and the velocity step size, Av, which is 
added in each compensation loop is 80 m/sec. The unknown target velocity, which is 
assumed between Mach 2 and Mach 4, is generated by using function "rand" in MATLAB. 
The simulation results are shown in Figure 3.26 through Figure 3.34. 

As determined from the simulation results of compensation scheme 3, the 


compensated HRR profile of the eighth loop (Figure 3.34) has the best target resolution 
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and the largest SNR. Thus, it is estimated that the relative radial velocity is approximately 


equal to 1220 m/sec. 
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Figure 3.26 The Uncompensated HRR Profile. The Relative Radial 
Velocity is Unknown, But between Mach 2 and Mach 4. 
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Figure 3.27 The Compensated HRR Profile of The First Loop. The 
Clutter is Canceled. Compensation Velocity = 660 m/sec. 
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Figure 3.28 The Compensated HRR Profile of the Second Loop. 
The Clutter is Canceled. Compensation Velocity = 740 m/sec. 
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Figure 3.29 The Compensated HRR Profile of the Third Loop. The 
Clutter is Canceled. Compensation Velocity = 820 m/sec. 
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Figure 3.30 The Compensated HRR Profile of the Fourth Loop. 
The Clutter is Canceled. Compensation Velocity = 900 m/sec. 
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Figure 3.31 The Compensated HRR Profile of the Fifth Loop. The 
Clutter is Canceled. Compensation Velocity = 980 m/sec. 
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Figure 3.32 The Compensated HRR Profile of the Sixth Loop. The 
Clutter is Canceled. Compensation Velocity = 1060 m/sec. 
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Figure 3.33 The Compensated HRR Profile of the Seventh Loop. 
The Clutter is Canceled. Compensation Velocity = 1140 m/sec. 
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Figure 3.34 The Compensated HRR Profile of the Eighth Loop. 
The Clutter is Canceled. Compensated Velocity = 1220 m/sec. 
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IV. CONCLUSION 


In the stepped-frequency radar, the high range resolution is achieved by performing 
the DFT on the received target signal. Nevertheless, this is optimum only when there is no 
relative radial velocity between the radar and the target. The nonlinear phase shifts 
induced by the relative motion cause a mismatch in the DFT process and result in a 
distorted range profile. Direct application of the velocity compensation factor to the 
received signal allows the compensation of the target signal; however, the clutter signal 
will disperse in the range profile due to the unnecessary compensation velocity. Thus, 
clutter cancellation has to be performed before velocity compensation in a velocity 
compensation scheme. In this thesis, three velocity compensation schemes have been 
presented to solve this problem. 

In the design of the velocity compensation schemes, weighting functions are used to 
suppress range sidelobes in the compensated range profiles. It is found that in all three 
schemes, the use of a Hamming window prior to the first DFT has the maximum influence 
on sidelobe reduction. In addition, the second compensation scheme was originally 
expected to have higher computational speed than the first scheme due to fewer 
processing stages. However, because of the circular convolution of the compensation 
factor and the target signal, the complex multiplies used in the second scheme are much 
more than those required in the first scheme, resulting in significantly longer special 
time. The first and the second schemes are theoretically equivalent, and both schemes 


generate very similar compensated HRR profiles for the same compensation velocity. 
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Nevertheless, since fewer Hamming windows are used in the second scheme, it is found 
that the compensated HRR profiles generated by the second scheme have higher range 
sidelobes and noise ripples. These phenomena can result in increased false alarms or 
reduced probability of target detection. From the viewpoint of fidelity and computational 
speed, the first scheme is considered better than the second scheme for implementation. 
The third scheme is based on the first scheme with the addition of a velocity compensation 
loop. In this case, the target velocity is unknown, but is assumed within a defined interval. 
In addition, it is observed that for the detection of moving targets using the 
stepped-frequency waveform, the use of a high PRF is appropriate for mapping high speed 
targets. The high PRF has the advantage of reducing the target dispersion even if accurate 
velocity compensation cannot be performed. However, a high PRF mode increases the 
system complexity. 

Waveform design for the stepped-frequency radar is much more complex than that 
for the constant frequency radar. The selection of waveform parameters involves several 
tradeoffs that are not present in standard radars. A study to highlight these tradeoffs and 


to streamline the waveform design process should be conducted in the future. 
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APPENDIX A. SIGNAL MODELING 


Target: 
In the scenario of the simulation, there are 4 dominant scatterers on the target. The 
magnitudes of N return pulses from a specific scatterer on the target are assumed constant. 


Thus, the sampled scatterer signal of the phase detector for the mth pulse is given as 


Ss; (n) =A; exp(/on) (A.1) 


where 


A, = the constant magnitude of a specific scatterer 
for the mth return pulse 


, = the phase component of a scatterer for the mth pulse 
n =1,2,3,.., N 


i = index/ is used to distinguished different 
scatterers on the target, i= 1, 2, 3, 4 


The phase component of the sampled scatterer signal for the mth pulse is given by Eq.2.7 


and Eq.2.9 as follows 
2Rn 
Dn _— 21 n Cc. 


~ ats ap(22== WWD) 
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4 — 1)A/)(Ro — (2 - 1)v2 
_ Amo +@ JAI o—@— Ivf) (A2) 


where 
Af = the frequency step size from pulse to pulse 


v = the relative radial velocity between the target and radar 
(- for opening velocity, + for closing velocity) 


T =pulse repetition frequency 
R, = the initial range when target is detected 


c =the speed of light, 3x10* m/sec 


Thus, a sampled scatterer signal is modeled as 


s; (n) =A; exp (Amo + WN Ro a DONRe _ dj) 


(A.3) 

Clutter: 

In the simulation, it is assumed that the stationary clutter distributes in an original 
range cell in which a moving target is located. Assume that there are 15 clutter sources, 
uniformly distributing in an original range cell (AR = 15 m) and the magnitudes of the 
return pulses from all 15 clutter source are constant. Similarly, the sampled clutter signal 


of each clutter source for the nth pulse can be modeled as 
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cx (n) = Cexp f a = aes) | (A.4) 


where 
C =the constant magnitude of the clutter signal 


R, = the range of each clutter source 


k =index kis used to distinguished different 
clutter sources, k = 1, 2, 3, ..., 15. 


Noise Source: 


The thermal noise is modeled as white Gaussian noise with variance o’. The 


signal-to-noise ratio (SNR) is given by 


A?/2 A? 


Thus, the constant magnitude of a scatterer signal can be written as 


A; = .[207(SNR) (A.6) 


The clutter-to-noise ratio (CNR) is given by 





C772 C* 
CNR = ae (A.7) 
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The constant magnitude of the clutter signal can also be written as 
C= 207(CNR) (A.8) 


Overall Received Signal: 
The overall received radar signal for the wth return pulse from an original range cell 


is modeled as 


S(”) = s S; (n) + > Cx (n) + noise(n) (A.9) 
i=] k=l 


where noise(n) is a white Gaussian sequence. 
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APPENDIX B. SIMULATION PROGRAMS 


function [ru, rstep]=schemel (N,fstep, prfr1,12,13,14,v) 
7/7/0700 000% %%%%%%0%0/070/0/0/0/0 
9%, | 

% This program is a MATLAB function which computes the compensated HRR profile of 
% 4-scatterer moving target in stationary clutter using scheme 1.This program include 

% clutter cancellation followed by velocity compensation. The input data consists of 

% following: 

% Number of pulses (N) 

% Frequency step size (fstep, Hz) 

% Pulse Repetition Frequency (prf, Hz) 

% Initial Scatter Ranges (r1,r2,13,r4, meters) 

% Target Radial Velocity (v, m/sec)(-opening)(+ closing) 

% The following parameters can be changed within the function 

% Signal-to-Noise Ratio in (dB) (SNR) 

% Clutter-to-Noise Ratio in (dB) (CNR) 

% Nominal Carrier Frequency (fo, Hz) 

% Velocity Matched for Compensation (co, m/sec) 

% The outputs are Ru=c/2*fstep, and Rstep=Kc/2*n*fstep, 

% which are the maximum unambiguous range and processed range bin size. 

% The FFT size is equal to N. 

% 

% 

TATA AV OOOO CUCL ALAA SASS SA MASA SA TATA SASA SATA SANA TATA SA TA 


% Radar Parameters 

fo=10*10%9; 

n=0:N-1; 

c=3*1048; 

f=fotn*fstep; 

K=1.33; 

PRI=1/prf; 

cnr=35; %dB 

snr=13; %dB 

a=10’(snr/10); 

A3=sqrt(2*a*0.0002); 

Al=sqrt(0.6)*A3; 

A2=sqrt(0.3)*A3; 

A4=sqrt(0.5)*A3; 

b=10(cnr/10); 

C=sqrt(2*b*0.0002); 
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% Signal Modelings, Target, Clutter and Noise 
s1=A1*exp(1j*(4*pi/c)*f.*(r1+n*v*PRI)); 
s2=A2* exp(1j*(4* pi/c)*f.*(r2+n*v*PRI)); 
s3=A3*exp(1j*(4*pi/c)*f.*(r3-+n*v*PRI)); 
s4=A4* exp(1j*(4*pi/c)*f.*(r4tn*v*PRI)); 
cl=C*exp(1j*(4*pi/c)* 1 *f); 
c2=C*exp(1j*(4*pi/c)*2*f); 
c3=C*exp(1j*(4* pi/c)*3 *f); 
c4=C*exp(1j*(4*pi/c)*4*£); 
c5=C*exp(1j*(4*pi/c)*5*f); 
c6=C*exp(1j*(4*pi/c)*6*f); 
c7=C*exp(1j*(4*pi/c)*7*f); 
c8=C*exp(1j*(4*pi/c)*8*f); 
c9=C*exp(1j*(4*pi/c)*9*f); 
c10=C*exp(1j*(4*pi/c)* 10*f); 
cl11=C*exp(1j*(4*pi/c)*11*0); 
c12=C*exp(1j*(4* pi/c)*12*f); 
c13=C*exp(1j*(4* pi/c)*13*f); 
c14=C*exp(1j*(4*pi/c)* 14*f); 
c15=C*exp(1j*(4*pi/c)* 15 *f); 
omegal=rand(size(n)); 
omega2=rand(size(n)); 
noise =sqrt(0.0002)*(omegal1+1j*omega2); | 
s=sl+s2+s3+s4+(cl+c2+c3+c4+c5+c6+c7+c8+c9+c10+cl1 1+¢12+c13+c14+c15)+noise; 
h=hamming(N)’; 
% Uncompensated HRR profile 
s=s.*h; 
S=ft(s); 
ru=c/(2*fstep); 
rstep=(K*c)/(2*N*fstep) 
figure(1) 
plot(n,20*log10(abs(S))) 
grid 
title((Uncompensated HRR profile ') 
xlabel(‘Relative range(meter)’) 
ylabel(‘Magnitude(dB)’') 
figure(2) 
t=linspace(1,ru,512); 
plot(t,20*log10(abs(S))) 
grid 
title((Uncompensated HRR profile ') 
xlabel(‘Relative range(meter)’) 
ylabel("Magnitude(dB)') 
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% Clutter Cancellation 

bu=input(Enter high bin number for clutter’) 
co=input('Enter compensation velocity for target ’) 
ve=exp(-1j*(4*pi*PRI/c)*(co)*f*n), 
ss1=[zeros(size(1:bu)) ones(size(1:N-bu))].*S; 
% Velocity Compensation 

SS 1=iffi(ss1).*ve; 

% Second Weighting 

SS1=SS1.*h; 

% Compensated HRR profile 

FSS=fft(SS 1); 

figure(3) 

plot(t,20*log10(abs(FSS))) 

grid 

title(‘Compensated HRR profile’) 
xlabel(‘Relative range(meter)’) 
ylabel(‘Magnitude(dB)') 

gtext('scatterers 1,2,3 and 4') 
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function [ru, rstep]=scheme2(N,fstep,prfr1,12,13,r4,v) 

OSS STA SSA NANA NANA NA SANA SATATANASATASANATANASASATASATATASALATA TAZA SA SATASA: 
% | 

% This program is a MATLAB function which computes the compensated HRR profile of 
% 4-scatterer moving target in stationary clutter using scheme 2.This program include 

% clutter cancellation followed by velocity compensation. The input data consists of 

% following: 

% Number of pulses (N) 

% Frequency step size (fstep) 

% Pulse Repetition Frequency (prf) 

% Initial Scatter Ranges (r1,r2,r3,r4) 

% Target Radial Velocity (v)(-opening)(+closing) 

% The following parameters can be changed within the function 

% Signal-to-Noise Ratio in (dB) (SNR) 

% Clutter-to-Noise Ratio in (dB) (CNR) 

% Nominal Carrier Frequency (fo) 

% Velocity Matched for Compensation (co) 

% The outputs are Ru=c/2*fstep, and Rstep=Kc/2*n*fstep, 

% which are the maximum unambiguous range and processed range bin size. 

“The FFT size is equal to N. 

% 

% 

AU SSSA SALSA SALAS SATA SA TASA SANA SANA SASANANASASATATASATASASANAZASATASATASASAS 


% Radar Parameters 
fo=10*1079; 
n=0:N-1; 
c=3* 1048; 
K=1.33; 
PRI=1/prf; 
cnr=35;%dB 
snr=13;%dB 
a=10’(snr/10); 
A3=sqrt(2*a*0.0002); 
Al=sart(0.6)*A3; 
A2=sqrt(0.3)*A3; 
A4=sart(0.5)*A3; 
b=10’(cnr/10); 
C=sqrt(2*b*0.0002); 
% Signal Modelings: Targets, Clutters & Noise 
f=fo+n*fstep; 
s1=A]1 *exp(1j*(4*pi/c)*f. *(rl+n*v*PRD); 
s2=A2*exp(1j*(4*pi/c)*f. *(r2+n*v*PRI)); 
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s3=A3*exp(1j*(4*pi/c)*f.*(3+n*v* PRI); 

s4=A4* exp(1j*(4* pi/c)*f. *(r44+n*v*PRI)); 

cl=C*exp(1j*(4*pi/c)*1*8); 

c2=C*exp(1j*(4*pi/c)*2*f); 

c3=C*exp(1j*(4*pi/c)*3*f); 

c4=C*exp(1j*(4*pi/c)*4*f); 

c5=C*exp(1j*(4*pi/c)*5*f); 

c6=C*exp(1j*(4* pi/c)*6*f); 

c7=C*exp(1j*(4*pi/c)*7*2); 

c8=C*exp(1j*(4*pi/c)*8*f); 

c9=C*exp(1j*(4*pi/c)*9*f); 

c10=C*exp(1j*(4*pi/c)* 10*f); 

c11=C*exp(1j*(4*pi/c)*11*f); 

c12=C*exp(1j*(4*pi/c)*12*f); 

c13=C*exp(1j*(4*pi/c)* 13 *f); 

c14=C*exp(1j*(4*pi/c)*14*f); 

c15=C*exp(1j*(4*pi/c)* 15 *f); 

omega 1=rand(size(n)); 

omega2=rand(size(n)); 

noise =sqrt(0.0002)*(omegal+1j*omega2); 

s=sl+s2+s3+s4+(clt+c2+c3+c4+c5+c6tc7+c8tc9+cl0tcl ]+c12+c13+c14+c15)+noise; 

h=hamming(N)’; | 

% Uncompensated HRR profile — 

s=s.*h; 

S=ffi(s); 

ru=c/(2*fstep); 

rstep=(K*c)/(2*N*fstep) 

figure(1) 

t=linspace(1,ru,N); 

plot(t,20*log10(abs(S))) 

grid 

title("Uncompensated HRR Profile’) 

xlabel(‘Relative range(meter)') 

ylabel(“Magnitude(dB)') 

figure(2) 

plot(t,20*log10(abs(S))) 

grid 

title(‘Uncompensated HRR Profile’) 

xlabel(‘Relative range(meter)') 

ylabel(‘Magnitude(dB)') 

% Clutter Cancellation 

bu=input(‘Enter high bin number for clutter’) 

co=input(‘Enter compensation velocity for target ’) 
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vc=exp(-1j*(4*pi*PRI/c)*co*f.*n); 
ss1=[zeros(size(1:bu)) ones(size(1:N-bu))].*S; 
% DFT of Velocity Compensation Factor 
ve=(1/N)*fft(vc); 

% Velocity Compensation 

SS 1=circonv(ss1,vc); 

% Compensated HRR profile 

figure(3) 

plot(t,20*log10(abs(SS1))) 

grid 

title((Compensated HRR profile’) 
xlabel('Relative range(meter)') 
ylabel(‘Magnitude(dB)') 

gtext('scatterers 1,2,3 and 4’) 
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function [value]=circonv(x1,x2) 
YH%%WHUHUHH%HH%%1%%%%7%%%%%U%%%%%%%%V0%/0%/0/0%0/0/0/0'/0 
% 

% This MATLAB function is to calculate cicular convolution via direct approach. Both x1 
% and x2 are sequences of length N. The resulting sequence of cicular 

% convolution has also the same length N. 

% 

% ) | 
AWWHRHUWUWUHUHVHUVUVUVLUVVUVUWUVVVUVWUWUUVVUWUVVV%%V%%V%0V0%0/0/07/0°/0/0 


k=2:length(x2); 
A(1,1)=x2(1); 
A(1,k)=x2(length(x2)-(k-2)); 
value(1)=sum(x1.*AC1,:)); 
for i=2:length(x2); 
for k=2:length(x2); 
A(i, 1)=AGi-1,length(x2)); 
AG,k)=A(i-1,k-1); 
end 
value(i)=sum(x1.*A(i,:)); 
end 
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function [ru, rstep]=scheme3(N,fstep,prf,r1,12,13,r4) 

LULA AoA >%% AVA AYA 
% 

% This program is a MATLAB function which computes the compensated HRR profile of 
% 4-scatterer moving target in stationary clutter using scheme 3.This program include 

% clutter cancellation followed by velocity compensation. The input data consists of 

% following: 

% Number of pulses (N) 

% Frequency step size (fstep) 

% Pulse Repetition Frequency (prf) 

% Initial Scatter Ranges (r1,12,1r3,14) 

% The following parameters can be changed within the function 

% Signal-to-Noise Ratio in (dB) (SNR) 

% Clutter-to-Noise Ratio in (dB) (CNR) 

% Nominal Carrier Frequency (fo) 

% Initial Compensation Velocity (v)(-opening)(+closing) 

% Velocity Matched for Compensation (co) 

% The outputs are Ru=c/2*fstep, and Rstep=Kce/2*n*fstep, 

% which are the maximum unambiguous range and processed range bin size. 

% The FFT size is equal to N. 

% 

% 

ASSL SASS STALLS LATA SA NASA SANA SASASANANASATATASATASASASATASASASALASANAS 


% Radar Parameters 
v=-(660+660*rand(1)); 
fo=10*10%9; 
n=0:N-1; 
c=3*10%8; 
K=1.33; 
PRI=1/prf; 
cnr=35; 
snr=13; 
a=10’(snr/10); 
A3=sqrt(2*a*0.0002); 
A1l=sqrt(0.6)*A3; 
A2=sqrt(0.3)*A3; 
A4=sqrt(0.5)*A3; 
b=10“(cnr/10); 
C=sqrt(2*b*0.0002); 
% Signal Modelings: Targets, Clutter & Noise 
f=fotn*fstep; 
s1=A1*exp(1j*(4*pi/c)*f. *(r1-n*v*PRD); 
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s2=A2*exp(1j*(4* pi/c)*f.*(12-n*v*PRI)); 
s3=A3*exp(1j*(4*pi/c)*f.*(r3-n*v*PRI)); 
s4=A2*exp(1j*(4*pi/c)*f.*(r4-n*v*PRD)); 
cl=C*exp(1j*(4*pi/c)* 1*£); 
c2=C*exp(1j*(4*pi/c)*2*f); 
c3=C*exp(1j*(4*pi/c)*3*£); 

=C*exp(1j*(4*pi/c)*4*f); 
c5=C*exp(1j*(4*pi/c)*5*f); 
c6=C*exp(1j*(4*pi/c)*6*£); 
c7=C*exp(1j*(4*pi/c)*7*f); 
c8=C*exp(1j*(4*pi/c)*8*f); 
c9=C*exp(1j*(4*pi/c)*9*f); 
c10=C*exp(1j*(4*pi/c)*10*£); 
c11=C*exp(1j*(4*pi/c)*11*f); 
c12=C*exp(1j*(4*pi/c)*12*f); 
c13=C*exp(1j*(4*pi/c)* 13*£); 
c14=C*exp(1j*(4* pi/c)*14*£); 
c15=C*exp(1j*(4*pi/c)*15*£); 
omegal=rand(size(n)); 
omega2=rand(size(n)); 
noise =sqrt(0.0002)*(omega1+1j*omega2); 
s=slts2+s3+s4+(cl+c2+c3+c4+c5+c6+c7+c8+c9+cl0+cll+cl2+c13+cl4tcl 5)+noise; 
h=hamming(N)’; 
% Uncompensated HRR profile 
s=s.*h; 
S=fit(s); 
mag=abs(S); 
ru=c/(2*fstep); 
rstep=(K*c)/(2*N*fstep) 
figure(1) 
plot(n,20*log10(mag)) 
grid 
title("Uncompensated HRR Profile’) 
xlabel(‘Relative range(meter)') 
ylabel(‘Magnitude(dB)’) 
figure(2) 
t=linspace(1,ru,512); 
plot(t,20*log10(mag)) 
grid 
title("(Uncompensated HRR Profile’) 
xlabel(‘Relative range(meter)’) 
ylabel(‘Magnitude(dB)') 
% Clutter Cancellation 
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bu=input(Enter high bin number for clutter’) 
co= -660; 

M=1:8; 

cco=co-(M-1)*80; 

for P=1:8 
vc=exp(1j*(4*pi*PRI/c)*cco(P)*f.*n); 
ss1=[zeros(size(1:bu)) ones(size(1:N-bu))].*S; 
% Velocity Compensation 
SS1=1fft(ss1).*ve; 

% Second Weighting 

SS1=SS1.*h; 

% Compensated HRR profile 

FSS=ffi(SS 1); 

figure(P+2) 

plot(t,20*log10(abs(FSS))) 

grid 

titleCompensated HRR profile’) 
xlabel(‘Relative range(meter)’) 
ylabel(‘Magnitude(dB)') 

end 
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